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We calculate nuclear suppression factor (Raa) f° r light hadrons by taking only the 
elastic processes and argue that in the measured px domain of RHIC, collisional rather 
than the radiative processes is the dominant mechanism for partonic energy loss. 

1. Introduction 

Jet quenching is one of the most promising tools to extract the initial parton density 
produced in high energy heavy ion collisions [H]. This is related to the final state energy 
loss of the leading partons causing depopulation of hadrons at high transverse momen- 
tum (jpt) ■ The suppressions of high px hadrons and unbalanced back-to-back azimuthal 
correlations of the dijet events in Au+Au collisions measured at Relativistic Heavy Ion 
Collider (RHIC) provide experimental evidence in support of the quenching. The ob- 
served nuclear suppression of light hadrons (tt, rj) in Au + Au collisions at ^Js = 62 — 200 
A GeV at RHIC was reproduced by radiative loss only assuming that the contributions 
from collisional loss is negligible. However, the non-photonic single electron spectrum 
from heavy meson decays measured by PHENIX Collaboration put this assumption 
in question. No realistic parameter set can explain this data using the radiative energy 
loss based jet tomography model which either requires violation of bulk entropy bounds 
or non-perturbatively large a s of the theory [ 4], or equivalently one requires excessive 
transport co-efficient q e ^ = 14 GeV 2 /fm 

In this light, we, in the present work would like to address if the omission of collisional 
loss at RHIC is justified or not. We argue that, whether the collisional or radiative 
loss is the main mechanism of energy dissipation, is a px dependent question. It also 
depends on the energies of the colliding system and expected to be different for RHIC 
and Large Hadron Collider (LHC). In contrast to the previous works, we calculate nuclear 
suppression factor for pions at RHIC energy considering only the collision energy loss. 

The partonic energy loss due to collisional processes was revisited in [E] and its impor- 
tance was demonstrated in the context of RHIC in[Cj. It is shown in ref.[E] that there 
exists an energy range where collisional loss is as important as or even greater than its 
radiative counter part, hence cannot be neglected in any realistic model of jet quenching. 
Recently this is also noted in ref.[HJ|B]. 
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2. Formalism 



We use Fokker Planck (FP) equation to dynamically evolve parton spectra. FP equation 
can be derived from Boltzmann equation if one of the partner of the binary collisions is 
in thermal equilibrium and the collisions are dominated by the small angle scattering 
involving soft momentum exchange [EHHIEIIE]- F° r a longitudinally expanding plasma, 
FP equation reads: 

where the second term on the left hand side arises due to expansion [H21- Bjorken 
hydrodynamical model [ IT5] has been used here for space time evolution. In Eq. (JTJ), /(p, t) 
represents the non-equilibrium distribution of the partons under study, r\ = (l/E)dE/dx, 
denotes drag coefficient, B\\ = d((Ap\\) 2 ) /dt, B± = <i((Ap_|_) 2 ) /dt, represent diffusion 
constants along parallel and perpendicular directions of the propagating partons. 

In our calculations we include some of the important features which were ignored in 
previous work [Ej in the context of jet quenching. These are: (i) the term correspond- 
ing to the longitudinal expansion is included, (ii) the momentum evolution of parton 
distributions both along longitudinal and transverse direction is considered and (hi) the 
mechanism of hadronization is introduced. 

The transport coefficients, rj, B\\ and B± appeared in eq.([TJ have been calculated in 
Ref. [[7j. The FP equation has been solved for the initial parton distributions parametrized 
as [[IS]: 

dN N o 

f(PT,P.,t = U) = j^y\y=° = TfTfp ( 2 ) 

where po = 1.75, a = 8. Since asymptotically the multiplicity of any hadron species 
produced by a gluon jet is 9/4 times that of quark [ ITBj . the normalization constant iVo 
for gluons is chosen accordingly. Solving the FP equation with the boundary conditions, 
f(p>t) — ► for \p\ — > oo, we evaluate the nuclear suppression factor, Raa defined as, 

„ / x J fa jp', T c )\ p > T=PT / z D a /^ (z, Q 2 )dz 

rCAA{PT> = „ , 7 / ; 71 n 7 n2U \P) 

IP ,Ti)\p' T =p T /z, JJa/nO{Z, y)dZ 

where /(p', Tj) and f(p',r c ) denote the parton distributions at proper time n and r c 
respectively. Here r, is the initial time and r c is the time when the system cools down 
to the transition temperature T c (=190 MeV). We have taken a s = 0.3 and the initial 
temperature Tj = 450 MeV. 



3. Results 

The result of our calculations for neutral pion is shown in Fig. ^ which describes the 
PHENIX data[HZ| for Au + Au at y/s = 200 GeV reasonably well. To stress our point 
further we also analyse the excitation function of the nuclear suppression factor. This 
has recently been studied in [[TH] attributing the suppression mechanism to radiative loss. 
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Figure 1. Nuclear suppression factor for Figure 2. Average parton energy versus 
pion. Experimental data are taken from transverse momentum of pion for ^/s = 200 
PHENIX collaboration [ 17 for Au + Au col- GeV/ A. 
lisions at ^/s = 200 GeV. Solid line indicates 
result from the present calculation with colli- 
sional energy loss of the partons propagating 
through the plasma before fragmenting into 
pions. 



As at lower y/s, collisional loss wins over its radiative counter part, we reanalyze the 
excitation function with the former and reproduce the data well [ITU]. 

To pin down the relative importance of 2 <-> 2 and 2 — > 3 processes, we determine 
the average energy of the parton which contribute to the measured px window of the 
hadrons. The average fractional momentum, (z) of the fragmenting partons carried by 
the pion is calculated using relevant parton distribution and fragmentation functions, i.e. 



(z) oc Jdzz{f(x a ,Q 2 )f(x b ,Q 2 )D a ^ h (z,Q 2 ) + f(x b ,Q 2 )f(x a ,Q 2 )D b ^ h (z,Q 2 )}. For the 



parton distributions, we use CTEQ [120] including shadowing via EKS98 parametrization[ 
l2*Tj . while for the fragmentation function KKP parametrization is used [1221- The average 
energy of the parton, (E parton ) is obtained by using the relation < E parton >= p T j < z > 
for y n = 0. The results are shown in Fig. |21 

It might be recalled that at RHIC energies the nuclear modification factor Raa{pt) has 
been measured in the pion transverse momenta range px ~ 1 — 13 GeV. Assuming that 
these pions are originated from the fragmenting partons it is clear that the maximum 
average parton energy required is about 26 GeV here. We might compare this value with 
the determined E c l given in ref. and note that at these energies collisional loss cannot 
be neglected. For lower beam energy, 62.4 (130) A GeV the value of maximum average 
parton energy required to produce a 13 GeV pion is 16 (22) GeV, where the collisional 
loss will definitely be more important. 

1 Note that E c is defined to be the energy below which elastic loss dominates [0. 
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4. Summary 

In conclusion, our investigations clearly suggest that in the measured pr range of light 
hadrons at RHIC collisional, rather than the radiative, is the dominant mechanism of jet 
quenching. Inclusion of three body elastic channels might even increase E c making our 
point stronger. 
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